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Abstract 

We introduce a new universality class of one-dimensional unimodal dissipative maps. 
The new family, from now on referred to as the (zi, Z2)-logarithmic map, corresponds 
to a generalization of the z-logistic map. The Feigenbaum-like constants of these 
maps are determined. It has been recently shown that the probability density of 
sums of iterates at the edge of chaos of the z-logistic map is numerically consistent 
with a (/-Gaussian, the distribution which, under appropriate constraints, optimizes 
the nonadditive entropy S q . We focus here on the presently generalized maps to 
check whether they constitute a new universality class with regard to (/-Gaussian 
attractor distributions. We also study the generalized (/-entropy production per 
unit time on the new unimodal dissipative maps, both for strong and weak chaotic 
cases. The (/-sensitivity indices are obtained as well. Our results are, like those for 
the z- logistic maps, numerically compatible with the q- generalization of a Pesin-like 
identity for ensemble averages. 
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1 Introduction 

One-dimensional nonlinear maps play an important role in the development 
of the theory of chaos. Their long-time behavior is different for different kinds 
of maps and, for dissipative dynamical systems, the phase space measure is 
not conserved: all trajectories approach a certain subset of the phase space 
called attractor. The characterization of chaotic attractors is interesting and, 
as we shall see, a still open problem whenever the Lyapunov exponent vanishes 
(frequently referred to as weak chaos). 

The exploration of their special dynamical properties is, besides their simplic- 
ity and convenience for the development of theory of chaos, also motivated 
by the hope that the study of the possible limits of validity of the canonical 
statistical mechanics can benefit from the study of much simpler dynamical 
systems that are known to exhibit statistical-mechanical analogies pQ. 

In particular, one-dimensional unimodal maps may depend of a single control 
parameter that determines the dynamical behavior of the map. They typically 
have only one attractor, which differs for different parameter values. This 
fact makes these maps to constitute paradigmatic models in the study of the 
emergence of complexity in dynamical systems. 

Here we introduce and analyze a new universality class of one-dimensional 
unimodal dissipative maps. Our initial scope is to test, through them, the ap- 
plicability and usefulness of generalized dynamical indicators (g-indices) that 
emerge within nonextensive statistical mechanics [2|3] in order to establish a 
more complete classification of (weak and strong) chaotic systems. 

The paper is organized as follows. In Sect.[2]we introduce the (zi, ^-logarithmic 
maps and we compare their attractors with those corresponding to two well 
known one-dimensional unimodal dissipative maps. In Sect.[3]we briefly review 
the generalized properties we are interested to test, as well as some peculiari- 
ties of their numerical study. In Sect. H]we present our numerical results. Our 
main conclusions are drawn in Sect. [5j 



2 One-dimensional unimodal dissipative maps: a new class. 

The well known z-logistic maps are among the simplest one-dimensional non- 
linear dynamical systems that allow a close investigation of complex behavior. 
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This family reads 

x t+1 = 1 - y\x t \ z (z>l;ne [0, 2]; \x t \ < V 



where /i is the control parameter, whose values are limited in order to avoid 
the orbits to escape to infinity; z characterizes the inflection of the map in 
the neighborhood of the extremal point x = 0. The larger is z, the flatter 
is the x = maximum. The z = 2 map is, as is well known, isomorphic to 
y t +i oc y t (l—y t ). These maps are known to have topological properties that do 
not dependent on z, and they constitute important universality classes of uni- 
modal maps. However, their metrical properties, such as Lyapunov exponents, 
chaos threshold control parameter values and Feigenbaum-like constants, do 
depend on z. In particular, the control parameter critical value \i c monotoni- 
cally increases from 1 to 2, when z increases from 1 to oo (we focus here, and 
hereafter, to the first entrance to chaos while increasing \x above zero). 

The ^-exponential family of maps was introduced j5] to characterize a fur- 
ther degree of flatness that z-logistic maps cannot attain [even for z — > oo in 
Eq. (pQ)]. They are inspired in Cauchy's exponential function (infinitely dif- 
ferentiable at x = and nevertheless nonanalytic). This family is defined as 
follows: 

xt+i = 1 - fie- 1 /' 9 *'" (z>0;iie[0,ii*(z)};\x t \<l), (2) 



where the upper limit fi*(z) of the control parameter guarantees that orbits 
do not diverge, and depends slowly from z [e.g., /x*(0.5) ~ 5.43]. It has been 
observed that there is a value of z above which the corresponding attractors 
are topologically isomorphic to those of the logistic map. These maps enabled 
to investigate chaos in a new universality class of maps [5]. However, their 
extreme flatness caused serious numerical problems when we tried to study 
various dynamical properties and, very particularly, the probability density of 
sums of iterates at the edge of chaos [5J. 

We then propose a new family of maps, to characterize chaotic behavior of an 
universality class different from that of unimodal maps. From now on we call 
them (zi, z 2 ) -logarithmic maps, and they are defined as follows: 



x t +i = 1 - H- 



x t 



[z x > 1; z 2 > 0; fi e [0, n*(z u z 2 )\] \x t \ < I) 



(3) 



where the parameters (zi, z 2 ) characterize the map. Notice that {z\, z 2 ) = (1, 0) 
{tent map) makes the map to be not differentiable at x — 0. The upper limit 
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fj>*(zi, Z2) varies slowly with (zi, z 2 ). These maps generalize the ^-logistic maps 
(which are recovered for z 2 = 0). At their extremum, they are less flat than 
the corresponding ^-logistic maps, in contrast with the ^-exponential maps. 

We have numerically verified, for a wide range of values of (z\, z 2 ), that the at- 
tractors of the (zi, ^-logarithmic map appear to be topologically isomorphic 
to those of the logistic map. As an example, Fig. [T] exhibits the /^-dependence 
of the dynamical attractor of the (1, l)-logarithmic map. This bifurcation dia- 




Fig. 1. Bifurcation diagram of the (l,l)-logarithmic map. 

gram shows a first whole cascade of period doubling of non chaotic attractors 
2 fc_1 — > 2 k at parameter values /i^, so that the sequence {nk\ (k = 1,2, • • •) 
approaches the critical value /x c = fi^ = 0.893425 . . ., the point that defines 
the chaos threshold. By further increasing fi above /x c , the behavior becomes 
complex and trajectories become chaotic with the exception of those 'win- 
dows' where stable periods re-emerge (see first window at // = 0.981 . . .). The 
entrances to such windows exhibit the intermittency phenomenon related to 
a tangent bifurcation, which implies the existence of a 'laminar phase'. Besides 
the period doubling scenario, the (zi, ^-logarithmic maps also exhibit band 
splitting phenomenon, so that when approaching \x c from the chaotic regime, 
a chaotic attractor band splits into two chaotic bands in such a way that 
the iterates alternate between both bands in a periodic way, even though the 
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Table 1 

Numerical values of the critical control parameter fi c ( z l, z "i)i an d of the Feigenbaum- 
like constant S(zi,z 2 )- In some cases, the period doubling sequence convergence is 
very slow. 
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6.19... 



movement is chaotic inside each band. A whole cascade of parameter values 
{fik} (k = 1,2, • • •) exists where there is a splitting from 2 fc_1 to 2 fc chaotic 
bands. Both the sequences of period doubling parameter values and the band 
splitting parameter values converge to the critical point \x c according to the 
equation: 

lim to ~ ^ = 5(z u z 2 ) , (4) 

fc->oo fi k+1 — fi k 

where 5{z\, z 2 ) is the Feigenbaum-like constant of the (zi, z 2 )-logarithmic map: 
see Tabled] We obtain analogous results for other (z%,z 2 ) values. Therefore we 
verify that, in spite of the fact that the topological properties do not depend 
on (zi,z 2 ), metrical properties do. The dependance of /i c on (zi, z 2 ) is depicted 
in Fig. El 
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Logistic map 



0.4 l 



Fig. 2. (z\,Z2) dependence of fi c . The red lines on the surface are also represented 
on a two-dimensional plot. The brown line represents the critical value /J, c (z) of the 
z-logistic map. 



3 Generalized properties of chaotic systems 



We are now ready to characterize chaotic attractors of the new family of 
maps. The z-logistic maps, and many others, have already been deeply studied 
[5ll8ll9lfT0] . However, some questions remain still open, especially at the edge 
of chaos. Furthermore, some already known properties might be tested on the 
presently introduced new family of nonlinear dynamical systems. 

In order to characterize chaos, we may say that a dynamical system is 'chaotic' 
if it possesses sensitive dependence on the initial conditions. That means that 
randomly chosen very close initial conditions generate totally different tra- 
jectories in the long-time run. When the temporal discrepancy increase is 
exponential, the system is said to be strongly chaotic (or just chaotic). But it 
is known that there are ubiquitous natural and artificial dynamical systems, 
typically at the edge of chaos, that exhibit a power-law sensitivity to the initial 
conditions. Such systems are termed weakly chaotic ones. In this case, stan- 
dard dynamical indicators give a poor description of the complexity of their 
time evolution. In particular, the Lyapunov exponent vanishes. Nonextensive 
statistical mechanics \3] provides a generalization of the standard dynamical 
indicators of chaos to characterize, not only strongly, but also weakly chaotic 
systems. 

Consistently, the first property we study is the sensitivity to initial conditions. 
The second property is the g-generalized entropy production per unit time (a 
concept that is devised to be analogous to the Kolmogorov-Sinai entropy rate 
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for strongly chaotic systems). In this work we are interested in studying these 
two properties on averages (denoted as (•••)) over the entire phase space of 
the maps (— 1 < x < 1) instead of using the quickest-spreading special regions 
[3]. This kind of averaging mimics Gibbs' approach to thermostatistical en- 
sembles, and minimizes the role of fluctuations (particularly large at the edge 
of chaos). The third property we are also interested in is the characterization 
of the probability density of sums of iterates of the map and, consequently, 
the possible applicability of a g-generalized Central Limit Theorem (CLT) to 
the iterates of deterministic dynamical systems. 

It has been shown in many one-dimensional unimodal maps that the sensitiv- 
ity to initial conditions is given by the generalized functional expression [3J 
(called g-exponential function within the context of nonextensive statistical 
mechanics) 

«')-Jj^„=-[ 1 + ( 1 -^V] 1/M - 3 V > <5) 



where Ax(t) is the temporal dependence of the discrepancy of two very close 
initial conditions at time t, q is a parameter (for q = 1 the expression recovers 
the standard exponential dependence e Alt ), and \ q is a generalized Lyapunov 
coefficient (when q = 1 then Ai = A and the Lyapunov exponent is recovered; 
when q < 1 and X q > the system is said to be weakly sensitive to the initial 
conditions; when q > 1 and X q < it is said to be weakly insensitive). 

We consider very close initial conditions, randomly chosen within the inter- 
val [—1,1], from which we estimate £(t). We perform this operation many 
times (typically 10 7 ), and average all the values of the corresponding \n q ^{t) 
(where ln q x = (x 1 ~ q — 1)/(1 — q) is the inverse function of the g-exponential; 
In! x = In x) for various values of q. We consider increasingly small initial 
discrepancies Aa;(0) between each pair of initial conditions such as to obtain 
results which no further depend on the value of Ax(0) for increasingly long 
times. We vary the value of q and verify a nontrivial property [5], namely 
that a special value of q, noted q™ n (where sen and av stand respectively for 
sensitivity and average), exists which yields a linear dependence with time. 
In other words, we verify that (ln 9 £)(t) ~ A"L t, where the linear coefficient 
Xgav constitutes a g-generalized Lyapunov coefficient. 

With respect to the entropy production per unit time, it is known that the 
Boltzmann-Gibbs entropy (Sbg = ~~ Y^Li Pi In Pi) is the appropriate one when 
strong chaos is present. But, at the edge of chaos (see [2.3.3] and references 
therein), we may conveniently use the g-generalized entropy 

i _ v w r> q 

s Z^iA {Sl = S BG ). (6) 
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For each universality class of maps (characterized by (z±, Z2)), one special value 
of the entropic index q, noted q ent (where ent stands for entropy), exists for 
which the entropy production is finite. We expect, at the light of many maps 
that have been previously studied, this entropic index to coincide with the 
one obtained from the study of the sensitivity to initial conditions (Eq. [SJ, 
i.e., q sen = q ent . Other methods (e.g., based on multifractality) do exist for the 
calculation of q, but we address here just these two ones. 

The estimation of the g-entropy production consists in dividing the phase 
space x in W (typically 10 5 ) equal cells, and putting N ic 3> W randomly cho- 
sen initial conditions inside one of the W cells. We accompany the spread of 
points within the phase space, and calculate S q (t) from the set of occupancy 
probabilities {pi(t)} (i = 1,2, ■■■,W). We repeat the operation many times 
(typically 10 3 — 10 4 for strong chaos, and W/2 for weak chaos), choosing differ- 
ent initial cells within which the N ic initial conditions are chosen (we usually 
use N ic = 10W). Finally, we average the entropies S q (t) over the N c initial 
cells so that the proper value of the entropic parameter q™ t is the special value 
of q which makes the averaged g-entropy (S q av) Nc production to be finite. The 
g-entropy production per unit time 

K q av = lim lim lim ^ S ^ Nc (7) 



is calculated taking into account that the partitions of phase space must be 
such as to obtain robust results. 

We also investigate the probability density of sums of iterates of the maps. The 
iterates of a deterministic dynamical system can never be completely indepen- 
dent, since they are generated by a deterministic algorithm. However, a Central 
Limit Theorem (CTL) for deterministic dynamic systems can be proved [7] 
when we consider a one-dimensional map, x t +\ = f{xt), with positive Lya- 
punov exponent. More precisely, the well known CLT assumption about the 
independence of N identically distributed random variables is replaced by a 
weaker property that essentially means asymptotic statistical independence for 
large time difference. In particular, the probability distribution of the rescaled 
sum 

1 N 

y = ^T,9(x t ) (8) 
jv t=i 



becomes a Gaussian for the number of iterates N — > 00, regarding the initial 
value random variable (7 = 1/2 for strongly chaotic maps). Here 

g : K d — > 3? fc is a suitable smooth function with vanishing average which 
projects from the <i-dimensional phase space to a fc-dimensional subspace. In 
our case, d — k — 1, and g(x t ) = x t — (x). It is rigorously proved that 
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the conditions of validity of a CLT, due to the mixing property associated 
with strong chaos, are satisfied for the logistic map, /i = 2. A CLT has not 
been rigorously proved neither for other parameter values nor for other z- 
logistic maps, but Gaussian limit behavior is also numerically observed in [H] 
for other z values on strong chaos regime. Consistently, we also expect to verify 
the Gaussian limit behavior on strong chaotic (z\, Z2) -logarithmic maps. It is 
clear, however, that this CTL does not hold at the critical points, where the 
Lyapunov exponent vanishes. Due to strong correlations between the iterates, 
a non-Gaussian limit behavior is expected in those points [9~fTU] . 

On the other hand, it is well known that g-Gaussian distributions 

p(y) cc e"^ 2 = (l + & (q - 1) y 2 ) (9) 

maximize the entropy S q under appropriate constraints. Consequently, on 
weakly chaotic z-logistic maps, where the g-generalized entropy (jSJ) is to be 
used, numerical indications of a g-generalized CLT are available [9ffTU] . With 
this respect, we are interested here in studying the distribution of the rescaled 
sums of iterates (|SJ) for the (z±, Zi ) -logarithmic maps. Notice that the rescaling 
factor iV 7 can be absorbed by calculating the variance a of the non-rescaled 
sum (7 = in Eq.fjHJ)) for a given N, and then plotting the histogram of the 
variable y/a. 

The sum (JSJ) must be evaluated for initial conditions located close to the 
space phase attractor. Consequently, we may omit the first iterates (i.e., a 
transient) until we obtain transient-independent distributions (a typical length 
of transient is up to 2 11 ). This is of course irrelevant in the N 00 limit, but 
it is numerically convenient when we must use finite values of N (typically up 
to 2 22 ). We use quadruple precision of Intel Fortran, to avoid roundoff induced 
effects. 



4 Numerical results 

We studied the sensitivity to initial conditions in both strong and weak chaos 
regimes. For strong chaos we obtain q a ^ en = q™ t = 1. This means that the 
Boltzmann-Gibbs entropy is the appropriate one for this regime. For weak 
chaos, instead, we obtain q™ n = q™ t < 1. 

These facts are illustrated in Fig.[3]for the (2, l)-logarithmic map. We obtained 
these entropic indices by fitting, over the intermediate regime (occurring before 
saturation), the curves with the polynomial A + Bt + Ct 2 and comparing their 
nonlinearity measure R = C/B. The optimum value of the entropic index 
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Table 2 

Numerical values, within the (z%, Z2 ) -logarithmic family of maps, of q™ n and q®% t 
indexes, on the Feigenbaum attractor (weak chaos). 

I „ „ \ „av „av 
\ z l, z 2) Isen lent 



(1,1) 
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± 





005 





36 ±0.01 


(1,2) 
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± 
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35 ±0.01 


(1,3) 
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± 





005 





34 ±0.01 


(2,0) 





358 


± 





005 





36 ±0.01 


(2,1) 





340 


± 





005 





34 ±0.01 


(2,1.5) 





336 


± 





005 





33 ±0.01 


(2,2) 





322 


± 





005 





31 ±0.01 


(2,2.5) 





315 


± 





005 





30 ±0.01 


(1-5,1) 





363 


± 





005 





38 ±0.01 


(2-5,1) 





321 


± 





005 





30±0.01 


(3,1) 





300 


± 





005 





31 ±0.01 



corresponds to R — (a straight line). The intermediate regime that we 
consider is such that the linear regression coefficient is constant (typically 
0.9999). Intrinsic fluctuations still persist, in spite of averaging (see Fig.3). 
We overcome them by studying the composed maps (j = 2,4). 

The same procedure is applied to other ^-logarithmic maps (see Table [2]). 
In all cases we obtain, within a small error bar, q™ n = q™ t . m Fig- HI we can 
see the influence of (zx, z 2 ) on the g-indices. 

Another interesting result that emerged is the coincidence, for both chaotic 
regimes, of the slopes of the sensitivity and entropy functions of time. For the 
(2, 1) -logarithmic map, we have iff = 0.382 ± 0.005 « Af = 0.372 ± 0.007 
(for strong chaos) and if^ 34 = 0.27 ± 0.01 « A^ 34 = 0.28 ± 0.01 (for weak 
chaos). These results reinforce those in [5], as they are numerically compatible 
with the g-generalized Pesin-like identity for ensemble averages. 

The probability distribution of the rescaled sums of iterates of the strongly 
chaotic (2, l)-logarithmic map presents, as expected, a Gaussian shape (see 
Fig. [S]). The weakly chaotic regime is more subtle. Indeed, it turns out to 
be necessary to gradually approach the exact value of \x c in order to attain 
the limit distribution. The critical parameter fi should approach /i c , while the 
number of iterates N of the sum (jSj) should diverge. In practice, N must be 
large enough to verify the N 00 TCL assumption, but not so large that the 
system "realizes" that \x is not exactly ji c . 
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Fig. 3. (2, l)-logarithmic map. Time dependence of the of (ln 9 £) (figures (a)-(b)) and 
(Sq) (figures (c)-(d)) for strong and weak chaos. The values of the control parameter 
are \i = 1.2 for strong chaos, and fi = 1.027082958927880 for weak chaos. Figure 
(e) shows the convergence of (Sq), for strong chaos, to the infinite fine grain limit 
(lent = !)• The (/-entropy function (S q ) averages over N c = W/2 cells. See that series 
of finer partitions in phase space do not change the value of the g-entropy production 
per unit time but, when increasing W, the first stage remains for a longer time: the 
rate of (/-entropy production for far-from-equilibrium evolution depends neither on 
the number ./V of points of the initial ensemble, nor on W. Figure (f) shows the 
convergence of (m 9 £), for weak chaos, to the infinitely small initial discrepancy 
of two different realizations (q™ n = 0.340). Sensitivity function (ln g £)(£) averages 
over Nr = 10 7 pairs of realizations; unless otherwise indicated, the discrepancy we 
consider is Ax(0) = 10~ 13 . 
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Fig. 4. (zi, ^-dependence of q™ en and q^t- These non-monotonic functions coincide 
within some small error bar. Dotted lines are guides to the eye. 



(2,1)-Logarithmic map 




-3-2-10123 

y/cj 

Fig. 5. Numerical estimation of the normalized probability distribution of the 
rescaled sums of iterates of (2, l)-logarithmic map, for two different control pa- 
rameter values (cr = 1.85 for fi = 1.05, and a = 6.62 for /i = 1.1). To study the 
problem near the attractor, we exclude a transient of Nt = 2 12 iterates. N c = 10 7 
is the number of random initial conditions and N = 2 12 is the number of iterates. 
Black line corresponds to the Gaussian function e - ^/' 7 ) , j3 = 3.1. 

First of all, we check that a transient time N t = 2 12 is enough to consider 
trajectories close to the phase space attractor. Omitting this transient, the 
distributions of the sums become independent of the transient length in all 
cases (see Fig. |6]). Summarizing, instead of El we use 

?/ = ^E g(xt) . (io) 

JV t=N t 
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Fig. 6. Numerical estimation of the normalized probability distribution of the 
rescaled sums of iterates of (2, l)-logarithmic map for the control parameter value 
u = 1.027083. We consider different transient times Nt- N is the finite number of 
iterates considered on the sum estimation. 



See, in Fig. [7J the probability density functions when we gradually approach to 
the exact value of \x c . An optimum intermediate value of N exists, for which 
the data collapse is produced. To obtain the 1.27-Gaussian convergence for 
an even more precise value of /x c , we expect that a much larger numerical 
value of N (N ^> 2 22 ) is needed. The numerical experiment becomes therefore 
computationally untractable. Fig. [8H9J illustrate the effect, on the distribution 
shape, of considering a finite value of N for two finite-precision values of 
the critical parameter \i c . In both cases, we find that intermediate values of 
N show g-Gaussian behavior. When N is not large enough (i.e., iV = 2 15 for 
H = 1.027083), the limit distribution exhibits a peaky shape in its central part, 
due to the fact that summation given by Eq. [8] is not adequate to approach 
the edge-of-chaos limiting distribution. 

Fig. [10] sketchily shows the effect, on the shape of the distribution, of consid- 
ering approximate values of the critical parameter with gradually improved 
precision, while iV increases. We then find a plateau of N values for which 
finite summation is adequate to approach the edge-of-chaos limiting distribu- 
tion, which plausibly is a g-Gaussian. 

This scenario is fully consistent with the one observed for the z = 2 logistic 
map [TP] . 
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Fig. 7. Data collapse of probability density functions for fi — ► [a c , and consistently 
enlarging N = 2 2n (see Fig. 10). 



(2,1)- logarithmic map 




-4 -2 2 4 



yP(0) 



Fig. 8. Rescaled probability density function at the edge of chaos with 4 digits 
precision (ju = 1.02709). For N = 2 18 it presents a 1.27-Gaussian behavior in its 
central part. For N » 2 19 , q — > 1, and probability density function tends to a 
Gaussian. 

5 Conclusions 



Let us summarize our main results: 

(i) A new universality class of one-dimensional unimodal dissipative maps is 
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Fig. 9. Approach to a q-Gaussian distribution with q = 1.27, for an approximated 
critical value with a 5-digit precision \i = 1.027083 (Au = 4 x 10" 8). 




( iH4 ) 



Fig. 10. Shape of the rescaled probability density function of a map at the 
edge of chaos, for finite N and finite precision for the critical parameter u c 
(u is the control parameter of the map). liniAr^oo lim^ Mc+ o PDF = peaky; 
lim (U _ >Atc +o hmAT-^oo PDF = Gaussian; limjv^oo PDF = q — Gaussian 
[fixed N(jjl — u c ) € (0, oo)]. By courtesy of the authors of Ref [10], where u = a, 
u c = a c , and s = log 4/ log 6. 
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introduced, characterized by a degree of flatness smaller than that of ^-logistic 
maps. 

(ii) The critical value on chaos threshold n c (zi,z 2 ) and the Feigenbaum-like 
constant are numerically estimated for (21,22) G [1,3] x [0,3]. Topological 
properties do not depend on (21, z 2 ) but metrical properties do. 

(iii) The entropic index q™ t , which makes the average of the g-entropy produc- 
tion finite, and q^ n coincide. This result is in accordance with the behavior of 
other classes of maps. For strong chaos we verify q™ t = q^ n = 1. 

(iv) The g-generalization of the Pesin-like identity is verified for ensemble 
averages for the (z±, 2 2 )-logarithmic maps, in both strongly and weakly chaotic 
cases. 

(v) Quantitative (zi, £ 2 )-dependence of the indices q™ t = q™ n is numerically 
studied for a certain range of (zi,z 2 ) values. 

(vi) The probability distribution of the sums of iterates in a strongly chaotic 
(zi, ^-logarithmic map is a Gaussian, as expected from the Central Limit 
Theorem for deterministic chaotic systems. 

(vii) The probability distribution of the sums of iterates for weakly chaotic 
(zi, ^-logarithmic map (i.e., at the edge of chaos) appears to approach a q- 
Gaussian, the probability distribution that maximizes the nonadditive entropy 
S q . These numerical results are consistent with a g-generalized Central Limit 
Theorem [TT] . 

All these results are expected to contribute to the correct interpretation of var- 
ious experimental features in dissipative dynamical complex systems [T2][T3"] . 
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